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I .  INTRODUCTION 


The  use  of  lasers  to  probe  the  combustion  process  has  become  a  very 
active  field  of  research.  Due  to  the  non-intrusive  nature  of  light  and 
the  monochromaticity  of  the  laser  this  tool  is  well  suited  to  make 
detailed  spatially  and  temporally  resolved  measurements  in  combustion 
environments.  Raman  scattering  techniques  (investigated  here)  can  be 
applicable  for  probing  the  combustion  process  at  elevated  pressures  since 
Raman  scattering  takes  place  in  less  than  10~1^  sec,  a  time  short  with 
respect  to  other  processes  such  as  collisions. 

We  report  here  the  use  of  spontaneous  Raman  techniques  to  measure 
the  flame  temperature  of  steady  state  flames.  The  CH4/N2O  flame  was 
investigated  most  extensively  in  preparation  for  studies  of  the  CH2O/ 

N2O  flame.  This  latter  flame  is  expected  to  be  important  in  the  flame 
zone  of  nitramine  propellants  since  N2O  and  CH2O  are  major  decomposition 
productsl,2  of  these  propellants. 


II.  EXPERIMENTS 

The  experimental  apparatus  used  for  generating  and  detecting  Raman 
signals  is  shown  schematically  on  Figure  1.  A  4  watt  cw  argon  ion  laser 
operating  at  488  nm  is  used  as  the  excitation  source.  This  laser  has 
its  cavity  extended  to  the  mirror  M  (2.5%  T,  5m  radius  mirror)  to  allow 
both  the  focussing  lenses  LI  (100  mm  focal  length,  AR  coated)  and  the 
sample  (slot  burner)  to  be  placed  inside  the  laser  cavity.  This  technique 
provides  a  significant  increase  in  laser  flux.  Ideally  one  could  obtain 
a  factor  of  ^  80  gain  using  the  2.5%  transmitting  mirror  for  the  laser 
output  coupler  M.  However  slight  losses  from  inserting  the  focussing 
lens  and  flame  inside  the  cavity  result  in  lower  gains.  The  laser  power 
is  monitored  by  measuring  the  portion  of  the  beam  that  is  transmitted 
through  M  with  a  thermopile.  The  thermopile  output  is  converted  to 
digital  data  and  recorded  with  a  single  counter  denoted  by  SC. 

Part  of  the  light  scattered  from  the  flame  is  collected  at  right 
angles  with  a  camera  lens  L2  (50  mm  focal  length,  fl.O)  and  is  focussed 
into  a  double  monochromator  DM  by  lens  L3  (254  mm  focal  length,  f5). 
Alignment  of  the  collected  light  with  the  monochromator  entrance  slit 
is  accomplished  with  a  dove  prism  D.  A  (thermoelectrically  cooled)  RCA 
31034  photomultiplier  is  used  as  the  detector. 
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Figure  1.  Schematic  Diagram  of  Raman  Apparatus 


Two  methods  of  recording  the  Raman  signal  are  used.  For  recording 
a  continuous  spectrum  the  output  from  the  discriminator  goes  through  a 
digital  to  analog  convertor  (DAC)  to  a  chart  recorder  where  the  time 
axis  generally  refers  to  scanning  the  wavelength  at  a  fixed  rate.  When 
background  emission  from  the  flame  is  high  another  method  which  subtracts 
the  signal  with  the  laser  off  (B)  from  the  signal  with  the  laser  on  (A) 
is  used.  This  is  accomplished  as  follows.  A  chopper  C  is  used  to  inter¬ 
rupt  the  laser  beam  and  a  trigger  signal  is  sent  to  a  clock  CL.  The 
clock  opens  the  gate  in  the  A  or  B  position  for  equal  amounts  of  time 
and  the  data  is  digitally  accumulated  with  a  dual  counter,  DC. 

A  microcomputer,  MC,  is  used  to  control  and  collect  data  when 
operating  in  the  gated  mode.  The  amount  of  time  desired  to  collect  data 
at  each  wavelength  together  with  the  wavelength  separation  between  data 
points  are  entered  into  the  microcomputer.  The  microcomputer  activates 
data  collection  for  this  period  of  time,  then  reads  the  data  from  the 
counters,  computes  A-B,  corrects  for  possible  laser  power  variations, 
and  computes  the  statistical  uncertainty  in  the  data.  These  values  are 
outputted  to  a  teletype  TTY.  Next  the  microcomputer  resets  the  counters 
and  changes  the  wavelength  of  the  DM  by  activating  the  stepper  motor  SM 
and  the  cycle  is  repeated. 

Flames  studied  in  this  work  were  produced  by  a  small,  water-cooled 
slot  burner  constructed  from  copper.  The  slot  dimensions  are  0.33  mm 
wide  and  28  mm  long.  Flow  of  the  fuel  and  oxidizer  to  the  burner  were 
individually  controlled  by  rotary  flowmeters.  In  the  results  reported 
for  the  slot  burner  measurements  were  taken  ^  2  mm  above  the  burner  head. 
The  spatial  resolution  of  the  sampled  signal  was  a  cylindrical  element 
65  ym  in  diameter  and  ^  1  mm  in  length. 

III.  THEORY 

The  Raman  effect  is  a  light  scattering  phenomenon  whereby  photons 
of  light  scattering  from  molecules  have  their  energy  changed.  If  the 
molecule  gains  energy  it  is  termed  Stokes  scattering  and  anti-Stokes 
scattering  if  the  molecule  loses  energy.  The  energy  changes  which  occur 
are  quantized  according  to  the  allowable  energy  levels  (rotational  and 
vibrational)  of  the  molecule. 

The  selection  rules  for  the  change  of  vibrational  (v)  and  rotational 
(J)  quantum  numbers  for  Raman  scattering  from  diatomic  molecules  are^>4 

Av  =  ±  1  and  AJ=0,±2.  (1) 
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If  one  visualizes  the  Raman  scattering  process  as  an  absorption-to  and 
emission-from  a  virtual  state  then  the  process  can  be  illustrated  as 
shown  on  Figure  2. 


A.  Rotational  Raman  Scattering 


Using  a  rigid  rotor  model  for  a  diatomic  molecule  the  rotational 
energy  eigen  values  are 


E 


rot 


h2-J  (J+l) 

8tt2  I 


(2) 


where  h  is  Planck* s  constant  and  I  is  the  moment  of  inertia  of  the 
molecule.  This  equation  is  usually  rewritten  as 


frot 

he 


=  F ( J)  =  BJ (J+l) 


(3) 


where  B  =  ^  is  called  the  rotational  constant  and  F(J)  is  called 

the  rotational  term  value.  Using  the  selection  rule  AJ=±2  the  frequency 
separations  (Av)  of  rotational  energy  levels  are 


Av  =  F ( J+2)  -  F (J)  =  4B (J+3/2) . 


(4) 


This  represents  a  series  of  Raman  rotational  lines  on  both  sides  of  the 
Rayleigh  line. 

A  rotational  temperature  can  be  calculated  from  determining  the 
rotational  quantum  number  J  which  has  maximum  intensity.  The  number  of 
molecules  in  rotational  level  J  of  a  ground  vibrational  state  at 
temperature  T  is 


N,  C2J+l)eBJ(J*1)hc/KT' 


J  occurs  when  — Hr  -  0. 
max  dJ 


which  gives 


(5) 
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Rayleigh  Rotational  Vibrational  Vibrational 

Stokes  &  Rotational  Rotational 

Anti-Stokes  Stokes  Anti-Stokes 

Figure  2.  Representative  allowed  scattering  processes  for  a  diatomic  molecule,  including  Rayleigh 
(elastic)  and  Raman  (inelastic)  processes. 


T  = 


(2J  + 

max 


1) 


2  hcB 
2K 


(6) 


For  rigid  diatomic  molecules  the  intensity  of  the  Raman  rotational 
lines  in  photons/s  is3-5 


R 


exp  [-BJi(Ji+l)  hc/KT], 


(7) 


where  Ji  and  Jf  refer  to  the  initial  and  final  rotational  quantum  numbers, 
Ar  is  a  collection  of  constants,  and  Qr  is  the  rotational  partition 
function  KT/hcB) .  The  light  frequency  of  Raman  scattered  light  is 
v  =  v0  ±  4B(J+3/2),  where  vG  is  the  frequency  of  the  excitation  source. 

The  ±  appears  in  order  to  include  both  Stokes  and  anti-Stokes  lines . 
S(Ji,Jf)  is  the  Raman  line  strength  which  is  given  as 


s(Ji,Ji+2) 


3 ( J i +1 ) (Jj+2) 

2(2Ji+3) 


(8) 


for  the  Stokes  case  and 


S(J.,J._2)  = 


3Jj(Jj-l) 

2(2J.-1) 


(9) 


for  the  anti-Stokes  case.  The  nuclear  spin  degeneracy  factor  gn  is  re¬ 
sponsible  for  the  even-numbered  J  rotational  lines  to  be  missing  for  O2; 
for  N2  the  even  J  levels  are  twice  as  intense  as  the  odd  levels. 

Rotational  temperatures  are  commonly  obtained  from  rearrangement  of 
equation  (7)  as 


F 

0 


(Ji)  ='  *n 


I(Jj,Jf)Ci(v) 

_S(Ji,Jf)gnv3  _ 


+  const. 


(10) 


** Michael  C.  Drake  and  Gera  M.  Rosenblatt,  " Rotational  Raman  Scattering 
from  Rremixed  and  Diffusion  Flames",  Combustion  and  Flame,  33_  179-196 


(1978). 
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and  plotting  such  that  hc/KT  is  the  slope  of  a  semi- log  plot.  Ci(v)  is 
the  calibrated  relative  intensity  response  for  the  experimental  detection 
system. 


The  effects  of  centrifugal  distortion  (i.e.,  a  non-rigid  rotor)  and 
contributions  due  to  rotational  Raman  scattering  from  vibrational ly 
excited  molecules  have  not  been  included. 


B.  Rotational -Vibrational  Raman  Scattering 

In  this  section  only  the  Q  branch  (AJ=0,Av=+l)  rotational-vibrational 
(ro-vib)  spectra  are  considered  for  a  diatomic  molecule  including  con¬ 
tributions  from  harmonic  and  anharmonic  oscillator  terms,  rigid  rotor 
terms,  and  ro-vib  interaction  terms.  The  individual  ro-vib  Raman  line 
positions  are  given  by 


V  =  Vo  ±  [T (v+1 , J)  -  T (v ,  J) ] . 


(ID 


H(v,J)  is  the  ro-vib  term  value  which  is 


+  [B-ag (v+1/2) ] J (J+l) . 


(12) 


We  is  the  equilibrium  vibrational  frequency,  WeXe  and  WeYe  are  the 
anharmonic  terms,  and  B  and  ae  are  rotational  constants.  Using  equation 
(12)  in  equation  (11)  we  have 


The  individual  line  intensities  are  given  by 


I(v,J)  a  gn(~2J+-Q')Q--+-^  v4C  exp  [  -  ^  (H(v,J))  ] 


(14) 


where  C  is  a  constant  and  0  is  the  vibrational  partition  function 


Q 

V 


1 
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For  a  comparison  of  experimental  data  with  equation  (14)  a  computer 
program  was  generated  to  produce  Q  branch  ro-vib  Stokes  Raman  spectra  for 
N2  (see  the  Appendix) . 


IV.  RESULTS 

One  of  the  first  gas  phase  pure  rotational  Raman  spectra  that  was 
recorded  was  obtained  using  a  cell  filled  with  100  torr  (1.33  x  104  pas¬ 
cals)  of  molecular  oxygen.  The  results  are  shown  on  Figure  3.  This 
spectrum  is  the  anti-Stokes  part  where  the  peaks  are  labelled  by  the  J 
rotational  state.  For  oxygen  only  odd  number  J  transitions  appear  due 
to  nuclear  spin  degeneracy.  The  intensity  of  this  rotational  spectrum 
is  maximum  around  J=9  from  which  an  equilibrium  rotational  temperature 
can  be  computed.  Using  equation  (6)  with  B  =  1.4456  cm-1  a  value  of  T  = 
375°K  is  obtained.  The  actual  room  temperature  of  the  cell  was  296°K. 
This  technique  is  not  very  precise  for  low  temperature  especially  using 
a  molecule  that  has  alternate  J  rotational  levels  missing.  However  it 
can  be  useful  for  obtaining  approximate  temperatures  with  a  minimum  of 
calculation. 

To  illustrate  another  method  of  obtaining  a  rotational  temperature 
a  rotational  anti-Stokes  Raman  spectrum  of  room  air  was  recorded  and  is 
shown  on  Figure  4.  Here  the  predominant  peaks  come  from  rotational 
transitions  in  02  and  N'2-  Only  the  rotational  levels  for  N2  separated 
from  the  02  levels  are  used  in  the  determination  of  the  temperature. 
These  levels  are  labelled  on  the  figure.  Using  these  intensities  in 
equation  (10)  a  temperature  is  obtained  from  the  slope  of  a  plot  of 


I(Jj,Jf)  C]  (v)  2(2Jj-l) 
3J.(J.-1)  gv3 


versus  Fo(Ji).  This  plot  is  shown  on  Figure  5.  An  unweighted  least 
squares  fit  to  the  data  gives  a  temperature  of  294  ±  15°K. 

This  method  has  been  used  to  measure  rotational  temperature  from 
flames  produced  on  a  slot  burner  operating  at  atmospheric  pressure.  An 
anti-Stokes  rotational  Raman  spectrum  of  a  rich  CH4-air  flame  was  re¬ 
corded  at  a  position  ^  2  mm  above  the  flame  reaction  zone.  Using  the 
N2  rotational  intensities  in  equation  (10)  a  rotational  temperature  of 
T  =  1230  ±  90° K  is  obtained  from  a  least  squares  fit  to  the  data  shown 
on  Figure  5. 

A  higher  temperature-flame  was  produced  by  using  a  near -stoichio¬ 
metric  CH4/N2O  mixture.  The  anti-Stokes  rotational  Raman  spectrum 
sampled  again  ^  2  mm  above  the  reaction  zone  is  displayed  on  Figure  6 
where  the  N2  peaks  used  to  calculate  the  rotational  temperature  are 
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Figure  3.  Rotational  anti-Stokes  Raman  spectrum  of  100  torr  of  oxygen  at  room  temperature. 
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Figure  4.  Rotational  anti-Stokes  Raman  spectrum  of  room  air 


Figure  5.  Boltzmann  plot  for  N2  from  rotational  Raman  spectra  recorded 
in  room  air  (+),  in  a  rich  methane-air  flame  (o),  and  in  a  near  stoichio 
metric  methane-nitrous  oxide  flame  (1). 
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Figure  6.  Rotational  anti-Stokes  Raman  spectrum  in  a  near  stoichiometric  methane-nitrous  oxide 
flame.  Nitrogen  transitions  used  in  temperature  calculations  are  identified. 
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labelled.  An  underlying  broad- smooth  peak  extending  out  to  200  cm~l 
is  believed  due  to  the  rotational  spectrum  of  CO2.  At  the  higher  rota¬ 
tional  J  values,  peaks  associated  with  the  rotational  spectrum  of  vibra¬ 
tional  ly  excited  N2  start  to  appear  out  from  under  the  ground  vibrational 
state  rotational  peaks.  The  intensities  obtained  from  this  spectrum  are 
again  cast  in  the  form  of  equation  (10)  and  plotted  on  Figure  5.  A  rota¬ 
tional  temperature  of  1830  ±  175°K  is  obtained  from  a  least  squares  fit. 

At  elevated  temperatures  excited  vibrational  states  become  suffi¬ 
ciently  populated  so  that  vibrational  temperatures  can  be  obtained  from 
a  comparison  of  the  relative  intensities.  The  cross  section  for  vibra¬ 
tional  Raman  scattering  is  not  as  large  as  that  of  rotational  scattering, 
but  it  is  usually  free  from  spectral  interference  of  other  molecular 
species.  Figures  7,  8,  and  9  show  Stokes  Q  branch  Raman  ro-vib  spectra 
for  N2  using  a  rich  CH4/air  flame,  a  near  stoichiometric  CH4/air  flame, 
and  a  near  stoichiometric  CH4/N2O  flame,  respectively.  The  peaks  are 
labelled  with  the  vibrational  quantum  numbers  of  the  states  participating 
and  the  asymmetric  broadening  of  the  peaks  is  due  to  rotational  transi¬ 
tions.  Vibrational  temperatures  are  obtained  from  computing  the  ratio 
of  the  relative  Q  branch  vibrational  intensities 


I  (v+1 , v) 
1(1  +  0) 


and  comparing  these  ratios  with  those  obtained  using  the  computer  program 
(Appendix)  with  the  appropriate  experimental  parameters.  The  vibrational 
temperatures  obtained  are  T  =  1180,  1680,  and  2180°K  for  Figures  7,  8, 
and  9,  respectively. 


V.  DISCUSSION 

Clearly,  at  low  temperature  pure  rotational  Raman  spectra  give  more 
accurate  temperature  measurements  since  excited  vibrational  state  popu¬ 
lations  are  small.  However,  as  the  temperature  increases  the  excited 
vibrational  population  grows  making  the  ro-vib  spectra  more  attractive. 

At  the  same  time,  rotational  transitions  within  vibrational  excited  states 
complicate  the  pure  rotational  spectra  resulting  in  a  decrease  in  the 
accuracy.  It  has  been  suggested^  that  below  temperatures  of  ^  2000 °K 
pure  rotational  Raman  is  more  accurate,  provided  rotational  spectra  of 
other  molecules  do  not  mask  the  spectra  of  interest.  Above  this  temper¬ 
ature  ro-vib  Raman  spectra  are  more  accurate. 


Q 

Michael  C.  Drake,  Chamnong  Asawaroenghai ,  and  Gerd  M.  Rosenblatt, 

" Temperature  from  Rotational  and  Vibrational  Raman  Scattering:  Effects 
of  Vibrational-Rotational  Interactions  and  other  Corrections ",  Div.  of 
Phys.  Chem.  ACS,  Washington  Meeting,  Sept  1979. 
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Figure  7.  Stokes  Q  branch  Raman  rotational -vi brational  spectrum  for  N2  in  a  rich  CIH4/air  flame  on 
small  slot  burner. 
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Figure  8.  Stokes  Q  branch  Raman  rotational-vibrational  spectrum  for  N2  in  a  near  stoichiometric 
CH4/air  flame  on  a  small  slot  burner. 
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Figure  9.  Stokes  Q  branch  Raman  rotational -vibrational  spectrum  for  N2  in  a  near  stoichiometric 
CH4/N2O  flame  on  a  small  slot  burner. 


Many  inves tigators5-15  have  utilized  Raman  spectroscopy  for  flame 
temperature  measurements;  temperatures  obtained  for  a  near  stoichiometric 
atmospheric  pressure  CH4/air  flame  have  ranged  from  1500  to  2000°K. 5 >7 ,9 , 15 
We  obtain  a  temperature  of  1680° K  from  ro-vib  Raman  spectra.  Adiabatic 
flame  temperature  calculations  from  a  version  of  the  NASA-Lewis  thermo¬ 
chemical  code!6  give  216915  and  2166°K  for  CH4/air  stoichiometric  mix¬ 
tures  at  atmospheric  pressure.  The  experimental  values  have  fallen  be¬ 
low  these  values  and  this  is  attributed  either  to  the  effects  of  heat 
loss  to  the  burner  head  or  to  cooling  from  the  entrainment  of  air  into 
the  flame. 

Both  pure  rotational  and  ro-vib  Raman  spectra  were  used  to  obtain 
temperatures  of  1830  and  2180°K,  respectively  for  near  stoichiometric 
mixtures  of  CH4/N2O  burning  at  atmospheric  pressure.  The  lower  temper¬ 
ature  obtained  by  pure  rotational  Raman  scattering  is  attributed  to  the 
influence  of  rotational  scattering  from  vibrationally  excited  molecules. 
Drake,  et  al.,3  calculate  that  this  effect  raises  the  temperature  by  8% 


7 

Alan  C.  Eckbrethj  " Laser  Raman  Gas  Thermometry AIAA  Paper  No .  74-1144 
Oct.  1974. 

g 

M.  Lapp ,  L.  M.  Goldman ,  and  C.  M.  Penney ,  "Raman  Scattering  from 
Flames ",  Science ,  175j  1112-1115  (1972). 

Q 

S.  M.  Schoenung  and  R.  E.  Mitchell ,  "Comparison  of  Raman  and  Thermo¬ 
couple  Temperature  Measurements  in  Flames ",  Combustion  and  Flame ,  35_ 
207-211  (1979). 

10 

W.  Strieker ,  "Local  Temperature  Measurements  in  Flames  by  Laser  Raman 
Spectroscopy ",  Combustion  and  Flame,  27,  133-136  (1976). 

11 

M.  Lapp  and  C.  M.  Penney,  eds..  Laser  Raman  Gas  Dia£nosticSj  (Plenum 
Press,  New  York,  1973). 

12 

D.  P.  Aeschliman,  J.  C.  Cummings,  and  R.  A.  Hill,  " Raman  Spectroscopic 
Study  of  a  Laminar  Hydrogen  Diffusion  Flame  in  Air" ,  J.  Quant.  Spectrosc. 
Radiat.  Transfer,  3_,  293-307  (1979). 

13 

William  M.  Arden,  Thomas  B.  Hirschfeld,  Stanley  M.  Plainer,  and  William 
A.  Mueller,  "Studies  of  Gaseous  Flame  Combustion  Products  by  Raman 
Spectroscopy ",  Appl.  Spectrosc.  28j  554-557  (1974). 

14 

A.  A.  Boiarski,  R.  H.  Barner,  and  J.  F.  Kircher,  "Flame  Measurements 
Utilizing  Raman  Scattering" ,  Combustion  and  Flame,  32,  111-114  (1978). 

15 

James  H.  Bechtel,  "Temperature  Measurements  of  the  Hydroxyl  Radical  and 
Molecular  Nitrogen  in  Premixed  Laminar  Flames  by  Laser  Techniques ", 
Applied  Optics  18,  2100  (1979). 

16 R.  A.  Svehla  and  B.  J.  McBride,  NASA  TND-7056  (1973). 
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for  N2  at  2000°K  using  a  2.6  cm-1  slit  width.  Our  slit  width  is  similar 
For  comparison  purposes,  the  temperature  has  also  been  measured  by  line 
reversal  using  the  OH  present  in  the  flame.  A  temperature  of  2211  ±  45K 
was  obtained.  A  3"  diameter  sintered  bronze  flat  flame  burner  has  also 
been  used  to  obtain  the  flame  temperature  for  a  near  stoichiometric  CII4/ 
air  atmospheric  pressure  flame. 1?  Initial  results  gave  temperatures 
around  2100°K;  more  recent  results  indicate  a  temperature  of  2200°K. 

Line  resersal  using  OH  was  also  used  with  this  burner  and  a  temperature 
of  2194  ±  25K  obtained.  The  calculated  adiabatic  flame  temperature  for 
this  flame  was  much  higher,  2922°K.18  Three  different  experimental 
techniques  suggest  the  flame  temperature  is  ^  2200° K.  Cooling  by  the 
burner  head  and/or  equilibrium  conditions  not  being  reached  may  account 
for  the  substantial  descrepancy  between  the  adiabatic  flame  temperature 

calculation  and  experimental  results. 


VI.  SUMMARY 

Extracting  a  temperature  from  both  pure  rotational  and  ro-vib  Q 
branch  Raman  spectra  of  molecular  nitrogen  in  a  flame  has  been  accom¬ 
plished.  A  flame  temperature  for  a  CH4/N2O  flame  has  been  measured. 
Error  limits  have  been  presented  for  rotational  Raman  temperature  data. 
Error  limits  for  ro-vib  Raman  temperature  data  will  be  addressed  else¬ 
where^  but  are  usually  less  than  ±  50K  for  clean  flames. 
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William  R.  Anderson >  " Laser  Excited  Fluorescence  Measurement  of  OH 
Rotational  Temperatures  in  a  CH4/N2O  Flame n ^  Eastern  Sectional  Meeting 
Combustion  Institute j  Atlanta 3  Georgia >  Nov  1979. 
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R.  A.  Hill,  A.  J.  Mulac,  D.  P.  Aeschliman ,  and  W.  L.  Flower,  "Temper¬ 
atures  from  Rotational-Vibrational  Raman  Q  Branches ",  J.  Quant, 
Speotroso.  Radiat.  Transfer,  21^  213-220  (1979). 

19 

R.  A.  Beyer  and  J.  A.  Vanderhoff,  " Raman  Spectroscopy  of  Premixed 
CH4/N2O  Flames ",  BRL  report  in  preparation. 
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APPENDIX 


This  Appendix  contains  the  computer  program  used  to  generate  the 
Stokes  Q  branch  ro-vib  Raman  spectrum  for  N2  as  a  function  of  temperature. 
Four  vibrational  levels  and  70  rotational  levels  are  included. 


MOLECULAR  AND  EXPERIMENTAL  PARAMETERS 


Symbol  in  Text 
W 


W  X 
e  e 


W  Y 
e  e 


a 


v 

J 


v 

o 


C 


Symbol  in  Program 


WE 

WEXE 

WEYE 

BE 

AE 

I 

J 

WO 

DEL  (slit  width 
triangular  function 
where  DEL  =  fwhm 

C 


Value 

2359.61  cm  1 
14.46  cm  1 
0.00751  cm 
1.9986 
0.01803 
0-4 
0-70 

20491.81  cm"1 

*  -1 
6  cm 

1X10"15 


This  program  calculates  the  0-1,  1-2,  2-3,  and  3-4  ro-vib  Raman  Q 
branch  spectrum  and  displays  the  plot  on  a  TV  monitor.  A  light  pen  is 
used  to  extract  the  peak  values. 


PROGRAM  LISTING 

C  DG  is  the  Raman  shift  and  WR  is  the  Raman 
C  frequency  for  the  ro-vib  Raman  line  positions. 
C  S  is  the  intensity  for  the  ro-vib  Raman 
C  Stokes  Q  branch  transitions.  ISUM  is  the 
C  intensity  convoluted  with  a  monochromator 

C  slit  function 

DIMENSION  WR  (6,80),  S(6,80),  ISUM(1500) 
EXTERNAL  PEN 

CALL  INTSET  ("324,  4,  ID,  PEN) 
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n  n 


CALL  DSPLPN  (1) 

REAL  I,  J 

C  ENTER  VARIABLES,  TEMPERATURE  AND 
THE  MONOCHROMATOR  SLIT  FUNCTION 
TYPE  100 

100  FORMAT  (3X,  'ENTER  T' ,  10X,  "ENTER  DEL'/) 

ACCEPT  11D,  T,  DEL 
110  FORMAT  (2F10.0) 

HCK  =  1.44 
WO  =  20491.81 
BE  =  1.9986 
WEXE  =  14.456 
WEYE  =  .00751 
AE  =  .01803 
WE  =  2359.61 
C  =  IE-15 
G  =  1.0 
K  =  1 

MIN  =  32000 
MAX  =  0 
H  +  DEL/6. 

Upper  and  lower  limits  on  frequency  for 
computation  of  synthetic  spectra. 

WL  =  WO  -  (WE  -  WEXE)  -  DEL 
WU  =  WO  -  (5* WEXE  +  DEL) 

NNX  =  4095/NPTS 
W  =  WL 

C  Rotational  and  vibrational  partition  functions 
QR  =  T/(2*HCK*BE) 

QV  =  1 ./ (1-EXP( (-HCK/T)* (WE  -  WEXE  +  .75  WEYE))) 

PRINT  120 

120  FORMAT  (3X  'WAVENUM' ,  13X,  'STOKES  Q  INTENSITY',  6X,  'V',  9X  'J'/) 

C  Calculation  of  the  ro-vib  Q  branch 
C  Stokes  Raman  line  positions  and  intensities 
DO  130  II  =  1,5 
DO  140  JJ  =  1,71 
I  =  II  -  1 
J  =  JJ  -  1 

DG  =  WE  -  2*WEXE*  (I  +  1)  +  WEYE*  (3*1*1*  +  6  +  13/4)  -  AE*J*(J+1) 

WR  (II,  JJ)  =  WO  -  DG 
F  =  INT  (J/2)  -  J/2 
IF  (F.LT.0)  G  =  .5 

S(II,  JJ)  =  (G/QR*QV)*  ( 2 * J  +  1)*(I  +  1)* (WR (II , JJ)**4) *C*EXP  (- (HCK/T) * 
(WE* (I  +  .5)  +  WEXE  * ( (I  +  . 5) **2)  +  WEYE*((I  +  .5)**3)  +  (BE  -AE/2)* 
J*(J+1)  -  AE*I*J* (J+l) ) ) 

PRINT  150,  WR(II, JJ) ,  S(II, JJ) ,  I,  J 
150  FORMAT  (IX,  E13.7,  10X,  E13.5,  7X,  F5.1,  5X,  F5.1) 

140  CONTINUE 
130  CONTINUE 
PRINT  160 
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160  FORMAT  (3X,  'MONWAVE'  13X  'SLIT  INT  INTENSITY'/) 

C  CONVOLUTION  OF  INTENSITIES  WITH 
MONOCHROMATOR  SLIT  FUNCTION 
170  CONTINUE 
SUM  =  0. 

SW  =  0. 

DO  180  II  =  1,5 
DO  190  JJ  =  1,75 
I  =  II  -  1 
J  =  JJ  -  1 

IF  (WR(II,JJ) .GE.W-DEL.AND.WR(II,JJ) . LE.W)  SW  =  (WR(II,JJ)  +  DEL  -  W) * 

S (II , JJ) /DEL 

IF(WR(II,JJ) .GT.W. AND. WR(II,JJ) .LE.W  +  DEL)  SW  =  (W  +  DEL  -  WR(II,JJ))* 
S (I I , JJ) /DEL 
SUM  =  SUM  +  SW 
190  CONTINUE 
180  CONTINUE 

PRINT  200,  W,  SUM 
200  FORMAT (IX, El 3 . 7 ,  10X,  EI3.5) 

IF  (W.GT.WU)GO  TO  230 
W  =  W  +  H 
I SUM  (K)  =  SUM 
K  =  K  +  1 
GO  TO  170 

C  SCALING  THE  DISPLAY 
DO  210,  M  =  1,  NPTS 
IF  (ISUM(M) . LT.MIN)  MIN  =  ISUM(M) 

IF  (ISUM(M) . GT.MAX)  MAX  =  ISUM(M) 

210  CONTINUE 

IDIFF  =  MAX  -  MIN 
IYSCALE  =  IDIFF/1023 
DO  220  M  =  1,  NPTS 
ISUM(M)  =  (ISUM(M)  -  MIN) /IYSCALE 
220  CONTINUE 

230  CALL  DSPMOl  (NPTS,  ISUM,  NNX,  0) 

GO  TO  230 
CALL  EXIT 
END 

SUBROUTINE  PEN  (ID) 

CALL  READ  (MOD,  IX,  IY) 

IF  (MOD. EQ. I)  TYPE  100,  IX,  IY 
FORMAT  (2 1 ID) 

IF  (MOD.EQ.l)  CALL  WAIT  (2.) 

END 
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TABLE  1.  RATIO  OF  0-1  TO  1-2  VIBRATIONAL  PEAKS  AS  A  FUNCTION 

OF  TEMPERATURE 


Temperature  DEL 

Kelvin  cm-1 


Vl-2 

VO-1 


1000  6 

1200  6 

1400  6 

1600  6 

1800  6 

2000  6 

2200  6 

2300  6 

2400  6 

2500  6 

2600  6 

2800  6 

3000  6 


0.078 

0.142 

0.218 

0.301 

0.389 

0.478 

0.565 

0.607 

0.647 

0.688 

0.728 

0.804 

0.816 
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01 


Temperature  versus  the  ratio  of  the  first  two  Raman  rotational - 
vibrational  peaks  of  N2,  from  computer-generated  Raman  spectra. 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


